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AN IMPROVED EXCEEDANCE THEORY FOR 


COMBINED RANDOM STRESSES 

By Harold C. Lester 
Langley Research Center 

SUMMARY 

This paper presents an extension of Rice's classic solution for the exceedances of 
a constant level by a single random process to its counterpart for an n-dimensional vector 
process. An interaction boundary, analogous to the constant level considered by Rice for 
the one-dimensional case, is assumed in the form of a hypersurface. The theory for the 
numbers of boundary exceedances is developed by using a joint statistical approach which 
fully accounts for all cross- correlation effects. An exact expression is derived for the 
n-dimensional exceedance density function, which is valid for an arbitrary interaction 
boundary. For application to biaxial states of combined random stress, the general theory 
is reduced to the two-dimensional case. An elliptical stress-interaction boundary is 
assumed and the exact expression for the density function is presented. The equations 
are expressed in a format which facilitates calculating the exceedances by numerically 
evaluating a line integral. The paper concludes with a brief discussion of the behavior of 
the density function for the two-dimensional case. 

INTRODUCTION 

In the last two decades considerable progress has been made in the application of 
random process theory to structural dynamics problems. Typical areas of application 
for some of these statistical methods have been (1) response of structural components 
such as shells and panels to boundary- layer turbulence and acoustical noise, (2) response 
of vehicles to surface (runway, rail, and road) roughness, and (3) response of aircraft and 
launch vehicles to atmospheric turbulence (wind gusts). With regard to this last area, 
that is, aircraft gust response, statistical methods have been developed to the extent of 
acceptance in aircraft design procedures. (See refs. 1 to 4.) Reference 1 contains an 
extensive bibliography. Contributions which are particularly noteworthy in this area are 
the advances which have been made in representing atmospheric turbulence as a continu- 
ous random process (ref. 5) and the accompanying refinement of power spectral methods 
for predicting aircraft gust response statistics (refs. 1 and 6, for example). Concurrently, 
advances in computer technology have made it possible to determine accurately and 



economically the significant flexible vibration modes of an airframe and to include these 
degrees of freedom in the structural model. 

In designing a structure, such as an airframe, to withstand the continuously varying 
stresses induced by a random loading environment, an important statistical response 
parameter is the number of exceedances. The exceedances are the expected (average) 
number of times per unit time that a single stationary random stress crosses a critical 
level (refs. 1 to 4, and 7). The solution to this classic problem in statistical mechanics 
was developed by S. O. Rice (ref. 8). Rice’s one -dimensional exceedance theory finds 
frequent application in fatigue -related studies, particularly for predicting service life 
under random loadings. The number of exceedances enables approximations to be made 
for the first-passage probability density function, mean time to failure, and the standard 
deviation of the time to failure. If the process is narrow band, then the probability dis- 
tribution for the peak stresses can also be calculated. (See ref. 7.) 

Another typical design situation (refs. 1 to 3, and 9) requires the consideration of 
two (or more) stress components acting simultaneously in combination. This situation is 
commonly referred to as a state of combined random stress. The individual components 
are interpreted as defining a random stress vector. For the two-dimensional case, the 
service life calculation is based on the exceedances, by the stress vector, of a closed 
curve which is called a stress-interaction boundary. 

Prior to reference 10, only the two-dimensional theory for combined random 
stresses had been developed. The first papers (refs. 2, 11, and 12) approximated the 
stress-interaction boundary as a connected union of straight-line segments. A question- 
able feature of these earliest exceedance theories is that they are based on relations 
which neglect some of the off-diagonal elements of the covariance matrix. The off- 
diagonal elements of this matrix account for all the cross -correlation effects between 
the stress and stress-rate components. Cooper (ref. 13) first recognized the potential 
importance of retaining the full covariance matrix and he accordingly derived an exact 
expression for the exceedance density function for a linear stress boundary. 

In the first part of this paper a general theory is presented for calculating the 
exceedances of interaction boundaries in an n-dimensional stress space. A joint- 
statistical approach, based on the full covariance matrix, is adopted. The stress bound- 
ary is taken to be a hypersurface. An exact expression is developed for the n-dimensional 
exceedance density function. This analysis is followed by the specialization of the general 
theory to two dimensions. The exact two-dimensional exceedance density is then pre- 
sented for an elliptical stress-interaction boundary. The equations for the density are 
expressed in a format which facilitates calculating the number of exceedances by numer- 
ically evaluating a line integral. (In ref. 10, an approximate closed-form expression is 
derived for this line integral by application of Laplace's method.) The paper concludes 


2 



with a brief discussion of the behavior of the density function for the two-dimensional case. 
Some useful properties relating to integration of quadratic forms are presented in the 
appendix. 

SYMBOLS 

ds element of surface area on stress -interaction boundary 

dX = dkjdxg- • -dx n 

dY - dy x dy 2 . . ,dy n 

E { } expected value 

Ej^ unit vectors, i = 1, 2, . . ., n 

erf(|) error function (eq. (40)) 


G(X),G(Y), 

G ( x 1> x 2)’ G < 0 ) 


exceedance density functions 


G(0) defined by equation (73) 

h(0) defined by equation (77) 

J Jacobian determinant 


K(Y) defined by equation (25) 

L,L orthogonal transformation matrices, equations (11) and (17), respectively 

element in ith row and jth column of L matrix 

M submatrix of M 


M covariance matrix (eq. (10)) 

M submatrix of Py (eq. (31)) 
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N average number of times per unit time stress vector crosses stress- 

interaction boundary 

N + exceedances, average number of times per unit time stress vector crosses 

stress -interaction boundary in outward sense 

n number of stress components 

P submatrix of M~ * (eq. (23b)) 

P submatrix of P~* (eq. (32)) 

p(X;X) joint probability density function 

Q submatrix of M~ * (eq. (23b)) 

Q submatrix of P~* (eq. (32)) 

R submatrix of M” 1 (eq. (23b)) 

R submatrix of P' 1 (eq. (32)) 

S(oj) power spectrum 

S submatrix of P y (eq. (31)) 

T submatrix of M (eq. (23a)) 

T submatrix of P y (eq. (31)) 

t time 

V vector (eq. (26)) 

V scalar (eq. (34)) 

V n component of stress-rate vector normal to stress-interaction boundary 

V| ith component of V y 
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w 

w n 

X(t) 

X(t) 

X'(t) 

x i 

*i(t) 

x lc’ x 2c 

Y(t) 

Y(t) 

yi 

y t (t) 

z 

z 

z i 

a(X),a(Y) 

a(d) 

a H 

*1 


vector defined by equation (A5) 
component of vector W 
stress vector, subvector of X(t) 
augmented stress vector (eq. (8)) 
vector (eq. (55)) 

coordinates of stress space where i = 1, 2, . . n 

ith component of X(t) where i = 1, 2, . . n 

critical stress levels for Xj(t) and X 2 (t), respectively 

stress vector, subvector of Y(t) 

augmented stress vector 

stress coordinates where i = 1, 2, . . ., n 

ith component of Y(t) where i = 1, 2, . . ., n 

vector (eq. (30)) 

subvector of Z (eq. (30)) 

ith component of Z where i = 1, 2, . . ., n 

defined by equation following equation (54) and equation (42), respectively 
defined by equation (74) 

element in ith row and jth column of P matrix 

root mean square of ith stress-rate component where i = 1, 2, . . ., n 
element in ith row and jth column of P" * matrix 



r 

correlation coefficient for x^(t) and Xgtt) 


V 

dummy variable of integration 


e 

coordinate defined by equation (69) 


m 

defined by equation (78) 


M 

mean stress vector (eq. (55)) 


Mi 

ith component of /I where i = 1, 2, . . n 


V 

defined by equation (76) 


P 

correlation coefficient for x^(t) and Xg(t) 


p' 

*0 

correlation parameter for x,(t) and x 0 (t) and also for x.,(t) and x 9 (t) 
(eq. (59g)} 

a. 

1 

root mean square of ith stress component where 

i = 1, 2, . . ., n 

d(8) 

defined by equation (75) 


*(U) 

defined by equation (A2) 


ct>(6) 

defined by equation (72) 


^/(X),^(x 1 

,Xg) stress-interaction boundary expressed as 

Mx) = o 

U) 

frequency, radians per unit time 



Capital letters generally denote vectors or matrices. Bars over Greek letters 
denote vectors. ( denotes matrix transpose. ( )"^ denotes matrix inverse. (’) 
denotes a differentiation with respect to time t. ( ) denotes ( ) referred to X space 

coordinates. ( ) denotes ( ) referred to Y space coordinates. ( ) + indicates that 

>y 

only outward crossings are counted. 

REVIEW OF EXCEEDANCE THEORY 

This section reviews previous one- and two-dimensional exceedance theories and 
notes two needed results from the literature. 
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Rice’s One-Dimensional Theory 


Consider, as shown in figure 1, a typical randomly varying function of time x^(t), 
for example, a stress induced in a structure by a random loading. It is a useful con- 
cept to think of this stress time history as a single realization of a large collection or 
ensemble of random stress time histories. Let x^ c be a constant critical stress 
level as illustrated by the dashed line. Of interest is the expected number of times 
per second N + that the random stress crosses the critical level in the upward sense as 
shown, for example, at points a and b. 


Xj(t) 



Figure 1.- Typical random stress time history showing exceedances 
of critical stress level x^ c - 

Rice (ref. 8) found that for a stationary normal random process with zero mean, 



in which x is the root- mean- square level of the stress x^(t) , and /3j x is the root- 
mean- square level of the stress rate x^(t). The root- mean- square values x and 

0^ x may be determined from the power spectral density S(a>) of the process x^(t) by 
integrations as follows: 


°1, X = J s ( w ) da) 


(2a) 


0? x = f w 2 S(w) dw 

“Loo 


(2b) 


The statistical average N + is commonly referred to as the exceedances and is a con- 
stant independent of time by virtue of the assumption of stationarity. 
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Two-Dimensional Theory 


When two or more random stress components act in combination, a state of com- 
bined random stress is said to occur and, in general, Rice’s one -dimensional exceedance 
theory is no longer applicable. An appreciation of the two-dimensional case may be 
acquired with the aid of figure 2. As shown in the figure, x^(t) and Xg(t) are two ran- 
dom stress time histories, perhaps a shear stress and an axial stress, respectively, and 
are interpreted to be the components of a two-dimensional stress vector X(t). When a 
design condition depends upon both stress components, it is possible to define a closed 

boundary curve \l/(x^,x^j =0 as illustrated in the figure. The closed boundary curve is 
referred to as a "stress-interaction boundary." In this case the exceedances N + are 
defined to be the expected number of times per second that the tip of the random stress 
vector crosses the boundary in the outward direction. Such crossings are illustrated by 
points a, b, and c of the figure. 



Figure 2.- Combined random stress showing exceedances of 
a stress -interaction boundary. 

Let N (note that the + superscript is missing) be the expected number of both 
inward and outward crossings. For the two-dimensional case, references 2, 10, 11, 
and 12 present expressions for this statistical average which are equivalent to 

N = ^ G^XpXg) ds (3a) 


/-k+OO ^1+00 

G ( x 1’ x 2) =)_ oo | V n|p( x l’ x 2* x l’ x 2) ^1^2 


(3b) 
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where 


(1) is the exceedance density function, crossings per unit time per unit 

of arc length on the boundary curve i//(xj,x 2 ^ = 0 


(2) <y ( ) ds indicates a line integral over the boundary curve 

(3) V n is the component of the stress-rate vector X(t) along the normal to the 

boundary curve 

(4) p( x i> x 2^i>^ 2) is i° int probability density function for the components of 

X(t) and X(t) 

EXCEEDANCE DENSITY FUNCTION IN n-DIMENSIONS 


The goals of this section are to note an extension of the preceding ideas to an 
n-dimensional stress space and to obtain an exact expression for the exceedance density 
function. 


The stress vector is assumed to have n components jxj(t), x 2 (t), . . x n (t)| i 
Cartesian space with a basis of orthogonal unit vectors |e^ x , Eg x , • . E n ,xj 


X(t) = Xl (t) E l x + x 2 (t) E 2 x 


+ x. 


En » J 


( 4 ) 


The components are random functions of time and are assumed to be jointly stationary 
with zero mean. Extension to random vectors with nonzero mean will be discussed 
subsequently. 

The arguments ordinarily given to support equations (3) may be paralleled by use of 
the geometry of n-dimensional vector spaces. In an n- space the stress-interaction 
boundary becomes a hypersurface 


if/(x v x 2 , . . x n ) = 0 (5) 

The exceedance density function and joint probability function take the forms 

G(X) = G^Xj, x 2 , . . ., x n ) 


9 



p(XjX) — p(xj, X£, . . Xj^j x^, X2, • • ■) x n ^ 


The expression, analogous to equations (3), for the total number of expected crossings 
may then be expressed in the form (ref. 10) 


N = § G(X) 


ds 


(6a) 


p+<x> r>+°Q r*+°° 

G(X)=j j ...] 

OQ *^-00 - OO 


p(X;X) dX 


(6b) 


where 


(1) dX = dxjdkg. • -dx n 

(2) ^) ( ) ds indicates an integration over the interaction boundary, now a 

hypersurface 

(3) The expression |v n | appearing in equation (6b) is given by (refs. 14 and 15) 


9* 

+ x 2 + 


1 9Xj 
1 9i p\ 

i 8x l I 


9Xr 


' dif/ 

, 9x 0 


+ X 


•n 


d\j/ 

% 


' dip 

i9x„ 


( 7 ) 


The stress and stress-rate components will be assumed to satisfy a joint normal 
(Gaussian) distribution. It is convenient to define a new vector X(t) (call it an aug- 
mented stress vector) by 



By virtue of the assumption of stationarity, it can be shown that X(t) has zero mean 
(ref. 16, p. 316). Since X(t) has been assumed to have zero mean, it then follows that 
X(t) has zero mean. The joint normal probability density then takes the special form 
(ref. 17) 
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(9) 


p(X;X) 


(27r)~ n 



exp 


- -i X T M' 1 X 


In equation (9) M is the covariance matrix for X(t) defined by 


M x - EjX(t) X T 

where E{} denotes expected value. By this definition covariance matrices are sym- 
metric and nonnegative definite. As discussed, for example, in references 16 and 18, a 
singular covariance matrix implies that one or more linear combinations of the random 
variables ^Xj, X 2 , . . ., x n ^ is zero with probability 1. In such cases the problem 

should be reposed in terms of a lesser number of random variables. Throughout this 
study M x will be assumed positive definite. The covariance matrix is therefore non- 
singular and the inverse operation in equation (9) is legitimate. 

Consider an arbitrary point on the stress-interaction boundary and visualize the unit 
outward normal at this point. Now define a new vector Y with components {Yj, Y 2 > • • • 
produced by the transformation 

Y = LX (11) 




The transformation is constructed so that the y^-axis is parallel to the unit outward nor- 
mal. The matrix L is a square orthogonal matrix of order n, and its elements are 
functions of position on the boundary. The transformation may be thought of as a rota- 
tion of axes. The orthogonality of the matrix L may be expressed as follows: 



The components of X also transform according to equation (11) so that 


( 12 ) 


Y = LX 


(13) 


The normal component of the stress- rate vector can now be written as 


Vn = Yi 


( 14 ) 


An augmented vector Y is defined as 



and it readily follows that 

Y = LX (16) 

in which L is the 2n by 2n matrix: 


L = 



(17) 


It follows directly from equations (12) and (17) that L is an orthogonal matrix 



(18) 


Substituting equations (8), (9), and (16) into equation (6b) yields an expression for 
the exceedance density function in terms of the coordinates of the Y' stress space. This 
expression is 


l mJ 


+0 °(27T)- n ' 


exp 


1 ~T 1~ 
- — Y M v Y 

2 y 


dY 


where 


dY = dyjdyg. . .dy n 


(19) 


( 20 ) 
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(21) 


T 

M y = LM X L 

In determining equation (19), use was made of the fact that dX = dY which follows from 
the orthogonality of L (ref. 19). 

For outward crossings y^ > 0 and therefore jy^ =y^. Modifying equation (19) 
accordingly gives for G + (Y), the exceedance density for outward crossings, 

r+ co r+ oo r +oo -n ~T~-1„ . 

G + (Y) = \ \ . . . \ — y, exp - ± Y M v Y dY (22] 

•-'-oo ^-oo )i ~ i 2 y 

\ M.J L 


The covariance matrix My and its inverse My may now be partitioned (see 
eqs. (A13); also ref. 17, p. 5): 


My. 


My ! T y 

T^iSy 


(23a) 


.-1 S 


(23b) 


where the various submatrices are all of order n by n. 

In order to place equation (22) into a form for exact integration, repeated use will 
be made of the ideas leading to equation (A17). With the use of this equation, equation (22) 
can be rewritten as 

p+°° r +cc r + °° i /• \ T i /• \ 

G +( Y> = K(Y) L . . . y, exp - |(y - V y ) P^Y - V y )J dy^y.,. . ,dy n 


13 


where 


exp - - Y T M" 1 Y 

L 2 y 


v y = " p y R y Y 


and (eq. (A14c)): 


p y - Sy - T W T v 


The components of V y are denoted as y) v 2 y> • • •» v n> yV 

Let Z be the vector with components z n } defined by 

Z=Y- Vy 

Equation (24) then becomes 

G+(Y > = K(Y) f " f ” . . . f ” (*! + v l>y )ej- \ Z^zl dz i dz 2- • • 

i,y L 


Partition once more and write 


z = T 


M JT 


A T I A 

T 1 JS 


ii, y 

“12, y • • ’ 

“ln,y 

21, y 

“22, y ' • • 

“2n,y 

nl,y 

“n2,y ’ • ' 

“ nn ,y 
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p" 1 = 

y 


q- 1 : r 


R T Ip" 1 


• 0 ! 


• fir 


As indicated Zj, M, and Q” are scalars (order 1). The vector Z is of order n - 1. 

A A A A .1 

The matrices T and R are of order 1 by n - 1. The matrices S and P are 

square matrices of order n - 1. 

By use once again of the results of the appendix, this partitioning allows equa- 
tion (29) to be written as follows: 

° +(Y) - k(y > c c • ■ • r (zi + vi > p ( 4i) exp [■ * (2 

I,y \ n i L 


v) P _1 (z - v) dz x dz 2 . . .dz n 


where the quantity V is a scalar given by 


V = - PR T z h 


The multidimensional integral on dzgdz-j. . .dz n is performed first and this 
integral is (eq. (A8)) 


p+oo „ + oo 

p+oo 

"|(z - v) T p- 1 (z - v) 

J ’ 

. \ exp 

V - OO - oo 

•' - OO 



in which (eq. (A14d)) 


P = P y M 


- i-i p y| 

m = -I — — 

“ii. y 
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Substituting equations (35) and (36) into equation (33) produces 


G + (Y) = (2tt) 


n-1 

2 




dz. 


(37) 


The remaining two integrations on z ± may now be performed and the results are 


r; 


z i exp r^ — dz i = “li.v exp - 


- \ 2a U,y 


i,y 


>y \ 2a 


ii, y/ 


(38) 


p+°° 


exp 


i,y 


2a 


dz, = - 


ii, y 


2 1 


1277a 


ii, y 


1 + erf 


h y 




(39) 


The error function erf(£) is defined by the integral 

erf(£) = 4= exp(-r] 2 ) dr? 

\frr J 0 


(40) 


and can be found in tabulated or series form in the literature (ref. 18). 

Substituting equations (38) and (39) into equation (37) gives for the exact integral 


n+1 


G + (Y) = 


in which 


(277) 2 ^ 

“u,y 



l 

a(Y) = 

V l,5 


P a ll 


exp (- a 2 (Y)) + \|ia(Y)|l + erf [a(Y)]j 


exp(- 1 Y T M~ 1 Y T j 


(41) 


(42) 


In determining equation (41), use has been made of equation (26) and the following relation 
which comes from equation (A14d): 
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1 


(43) 



Kl N 

Equation (41) constitutes an exact expression for the exceedance density function 
and is valid for an arbitrary stress -interaction boundary. The density G + (Y) is a point 
function of the Y space stress coordinates. It is desirable, therefore, to reexpress the 
density function in terms of the original X stress coordinates. 

T - 1 

By substituting equation (11) into the homogeneous quadratic form Y My Y 
(eq. (41)), it is quite easy to show that 

Y T M^ 1 Y = X T M^ 1 X (44a) 




(44b) 


An alternate form of equation (26), which defines the vector Vy, will be used: 

V y = tJm^Y (45) 

This form follows from equation (A14a). The IV^ and M~* matrices (eq. (10)) are now 

^ ~ _ 1 

partitioned in a like manner to the partitioning of the My and My matrices in 
equations (23): 





nT! 



(46a) 



Qx 1 

Rx 

** 

p ; 1 


(46b) 
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Next, by using equations (10), (16), and (46), it can be shown that 


M y = LM X L 


T y = ltJl t 




Substituting equations (16) and (47) into equation (45) gives 


V v = LV X 


(47) 


(48) 


where 


V x = tJm x X X (49) 

The components of the vector V x are denoted as £v^ x , Vg x , . . v n x}- 

Let L j be the scalar element in the ith row and jth column of the transformation 
matrix L. The elements fjj (j = 1, 2, . . ., n) are the components of the unit outward 
normal vector in the X system. From equation (48) the first component of V y is 


n 


!,y 


•l 

i=i 


Vj,x 


(50) 


The remaining quantity to be transformed is 
the matrix P x are denoted by 


«ll,y (eq- ( 31 ))- 


The elements of 


P x ^x T x M x T x 


11, X 

“12, x • • 

• “ln,x 

21, x 

“22, x • • 

• “2n,x 

nl,x 

“n2,x ' • 

“nn,x 


(51) 


18 


I 



and, as before, 


Py = LP X L J 


(52) 


Substituting equations (31) and (51) into equation (52) gives 


n n 


“ 11 , y = 2 2 { 11 Vij,; 


i»l i=l 


By collecting results 


n+1 


G + (X) = 


(2tt) \a 


11, y 


M v 


exp^-Q!^(X)^ + ^7TQf(X){l + erf(a(X))} 


exp-lxX'x 


in which (from eqs. (42), (50), and (53)) 


(53) 


(54) 


iZ. 


a(X) = - 

f“n,y 

T l,y “ 2 Vi,, 

i=l 


n n 

“li, y = 2 2 J li‘ij“ij,x 
i=l j=l 

This section is concluded by commenting on the modifications when the random vec- 
tor process has nonzero mean. Assume now that X(t) has a nonzero- mean vector p x 
with components x , M 2 x’ ' ‘ •> ^n x}- The components of the mean vector are con- 
stants for stationary processes. Define a new random vector process 


X’(t) = X(t) - Ji 


(55) 
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x 2> E 2 ,x 



Figure 3.- Stress-interaction boundary showing effect 
of nonzero-mean stress vector. 


The process X'(t) has mean zero by definition. As illustrated in figure 3, for the two- 
dimensional case, the transformation equation (55) may be viewed as a shift of origin. 
Therefore, to apply equation (54) to this situation, it is necessary only to 

(1) Replace X(t) by X'(t) 

(2) Replace the covariance matrix M x by the covariance matrix M x , for the 

(zero-mean) X'(t) process 

(3) Express the stress-interaction boundary in terms of the x! stress coordinates 

as follows: 

^(x’ + n ljX , x' + (i 2>x , 4n,x) = 0 

THE TWO-DIMENSIONAL CASE 

In this section an elliptical stress-interaction boundary is assumed and the previous 
results are specialized to two dimensions. The characteristic behavior of the exceedance 
density function is then discussed. 

Reduction of n-Dimensional Theory 

For the two-dimensional case the augmented stress vector has the form 
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M 


X(t) = < 


x iW 

xj(t) 


x 2 (t) 


(56) 


in which Xj(t) and x 2 (t) have zero mean. As previously discussed, the bias imposed 
by nonzero-mean stress components is accounted for by expressing the stress-interaction 
boundary as 


<//(xi,x 2 ) = 


X 1 + p l,x\ x 2 + p 2 


^ + 


£ - 1 = 0 


“■lc 


4 2c 


(57) 


The constants x are the mean stresses for the original random processes 
x i(t) + M i;X (i = 1, 2)/ 

The covariance matrix for X(t) will be written as 


M x = 


a l,x pa l,x°2,x 

pa l,x a 2,x °2,x 


-pa, o„ 
'o l,x 2,x 


P Oi On 

H o l,x 2,x 
0 


O C 7, On 

l,x 2,x 


l»x 


-P CT, On 

H o l,x 2, X 
0 


y ^l,x^2,x 
2 


y Pl,x@2,x @2,x 


(58) 


olx = E x2(t) 


°2,x = E ( x 2^ 


Pl,x = 


:,2/ 


(59a) '• 
09b) 


(59c) 

i“ It. 


2T>' 


in which 



(59d) 



y = 


h,x^2,x 


x^t) x 2 (t)> 


(59e) 


(59f) 


p o = 


CT l,x a 2,x 




(59g) 


1 1 ^ 

The zero off-diagonal terms appearing in M x follow from the well-known results for 

stationary processes: 



The correlation coefficients p and y and the cross- correlation parameter p' Q can be 
shown to satisfy the following inequalities (ref. 16): 


-llpil 


-1 S y S 1 





(60) 


It should be noted that p' is the cross- correlation parameter which is neglected in the 

r; ^ f 

exceedance theories of references 2, 11, and 12. If p Q = 0, then the covariance matrix 
M x (eq. (58)) assumes a block- diagonal form and the calculation of the exceedances is 
considerably simplified. 

The elements St ^ (j = 1, 2) of the transformation matrix L are the components 
of the unit outward normal vector to the elliptical boundary: 
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x l + ^l,x 



x 2 + ^2,x 



(61) 


By use of equations (46a) and (58), 

f ,2 


1 


1-P 2 


la 


l,x 

■ p Al,x a 2,> 


P/ a l,x°2,> 


1 /cl 


2 ,x 




5 

J 

The determinant for M x is 



M x 

= 0 , cr„ (l - 
l,x 2,x\ K / 


From equations (46) and equation (A14c) 




<*?,x - 

vh,xh,x - pii \x°2,x 


p x = 

lh,xh,x - ri\x a 2,x 

4,x - ^.x . 


(62) 


(63) 


(64) 
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where 


f.SL 

1 - P 5 


Now from equations (46a), (49), and (62) 


. .JSlX 

1,x 1 - p 2 \ a 2,x , 


p o / CT 2,x ^ ^ 


Equation (50) becomes 


v l,y = * ll T l,x * { 12 v 2,x 


Next, with the aid of equations (51) to (53) and (64) 


a tl,y " (*l,x " + 2 (^l,x%,x " p ^ 2<7 l,x a 2,x) f 11^12 + ^2,x ' p2a 2,x| C 12 (68) 


In order to facilitate numerical integration a point on the stress-interaction bound- 
ary is expressed in terms of a single parameter 6 as follows: 


x l + p l,x = x lc COS 6 


x 2 + M 2;X = x 2c Sin 6 


The exceedances N + (see eq. (3a)) can now be determined from 


+ C 2ir + 
n + = \ G + (e) 

Pn 
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whereby using the transformation equations (69) and previous results, G + (0) can be com- 
puted from the following relations: 


a(6) = 


d{6) = 


h(0) = 


G + (9) = <p(6) exp(i/h(0)) 


(71) 


x lc V x 2c 




\ a l,xJ\2,xJ z 


(2,) 3/2 fT- 


G(0) o(6) 


(72) 


G(0) = exp(-a2(0)) + fjj a(d) 


1 + erf { 01 ( 0 )} 


(73) 


l/2(l - p 2 )d(6) 


- p cos 20 


[ 1 x 2c/ g 2,x 


2 x 


lc 




L lc 


£12 


4 2c 


/° 2 ,> 


sin 20 + 


/> 1. 


L lc 


- ^ ^/°2,xV 0S , . Lh±_ JUafis/^L 4 

x 2c X lc/ G l,x/ \ x 2c x lc x 2c/°2,xj j 


A,x/ CT l,x 


2 -p2 y ^£^£-pp2 

’ vy cos 2 0 + 2 — — sin 0 cos 0 + 


(74) 
1 1/2 




sin^ 0 


/x-i /CT ^ 

x lc x 2c / 

t lc / l,xj 

<*1 - ^ ' 


( x 2c/°2,: 


v = 1 

2(l - P 2 ) 

. \2 xx /x \ 2 

— — | cos20 - 2p — — — — sin 0 cos 0 + ^ — | sin 2 0 + 2(1 - P 2 )x(0) 

Or, / ff, a„ a o \ / ' 

v l,x/ l,x 2,x \ 2,x/ 


(75) 

(76) 


(77) 
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X(0) = 


1 - P* 


r 

x lc / m 2,x x 2c _ M l,x x lc 
CT l,x\ x 2c ct 2,x x 1c a l. 


v. 


2r 

cos 0 + Ip 

a 2 ,x 


M l,x x lc 


m 2,x x 2c 


x, a, 
lc l,x 


x 2c °2,x/ 


sin 0 


4 ' 



2 P 


M l,x m 2,x x 1c x 2c 


x lc x 2c CT l,x a 2,x 






( 78 ) 


Characteristic Behavior of Exceedance Density Function 

This section presents some trends in the behavior of the exceedance density func- 
tion for the two-dimensional case. The elliptical stress-interaction boundary will be 
specialized to a unit circle. The cross-correlation parameters p, y, and p' Q and the 
root- mean- square parameters a, a„ /3, . and j3 0 „ will be varied from the fol- 

X y di. y XV X y xV jX 

lowing base values: 


CT l,x = ct 2,x = 1 

01 ,x = 02, x = 1 
P = y = Pq = 0 \ 

^l,x = ^2,x = 0 
x lc = x 2c = 1 


(79) 


The main purpose is to indicate the manner in which these parameters individually affect 
the behavior of the density function. For this reason results will be shown for some 
parameter combinations which are not physically realizable. For example, with station- 
ary processes p = 1 implies that y = 1 and also that p' Q = 0, whereas results will be 
presented for p — 1, y = 0 and p = 0, y — 1. The bounds imposed by equations (60) 
will be maintained. 

Figure 4 shows the effect of varying p, the cross -correlation coefficient for the 
stress components x^(t) and Xg(t). For a value of p = 0, the density is a constant 
function. This condition is expected because with zero cross correlation and a circular 
boundary, there is no preferred direction for the stress vector. As p approaches unity 
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(P - 1)> prominent peaks emerge in the density at 9 = 45° and 8 = 225°. This is also 
expected. As p - 1 the stress vector tends to fall more and more along the line 
x 2 = X 1 s h° wn i n the inset, top right of figure 4. For the limiting value of p = 1, it caj 
be shown (ref. 16) that ^(t) = Cxj(t) with probability 1. For the present case, 
the constant C = 1. It thus appears plausible that as p - 1, the density would take on 
the characteristics of two delta functions located 180° apart. 



Figure 4.- Exceedance density function showing effect of cross- correlation 
coefficient p. (For base case a 1>x = * 2 ,x = h,x = h,x = U Y = p' 0 = 0; 

Pl,x = H,x = °5 x lc =x 2c = 1 -) 


The inset plot, top center of figure 4, shows the variation of the number of exceed- 
ances N with p. The numerical value of N + was obtained by numerically evaluating 
the area under the corresponding density curve. For values of p up to about 0.6, the 
curve is flat and then gradually falls off for higher values of p. 

Figure 5 shows the bias imposed on the density by unequal root- mean- square stress 
levels. Here an off-nominal correlation coefficient of p = 0.6 has been arbitrarily 
selected and the densities plotted for a. =1 and ^ x = 2. As might be expected the 
peak at 45 for Oj x = ^ shifts toward tne x^-axis when o ^ = 2. A similar shift is 

noticed for the peak originally occurring at 225°. Exceedance values for the two cases 
are shown in the table at the top of the figure. Doubling a « while holding the root- 



mean-square stress rate f}-. fixed has the effect of reducing the number of exceed- 

i »x 

ances. The same trend is apparent in Rice’s one-dimensional theory. (See eq. (1).) 



Figure 5.- Exceedance density function for p = 0.6 showing effect of 
root-mean-square level Uj x - (For base case, a 2 x = ^1 x = 


y 

o 


^2,X " ^ - ^1 y ^ 9 . y ~ 


l,x ^2,x 


x lc = x 2c 


-*•) 


In figure 6 the effect of varying y, the cross-correlation coefficient for the stress- 
rate components x^(t) and ^(t), is shown. As the correlation is increased from y = 0 

(the base case for which G + (0) is a constant function) the density again begins to peak at 
45° and 225°. The peaks, however, are not as sharp as in the previous cases (fig. 4). 

This behavior is also expected. For y = 1, reflecting perfect correlation between the 
stress-rate components, X 2 = x^ with probability 1. Therefore, the velocity vector 
remains parallel to the line X 2 = Xj. For this limiting situation the density is zero only 
at the tangent points of 135° and 315°. 

The exceedances N + for this case are shown plotted against y in the inset of 
figure 6. The behavior is similar to that discussed previously. 

The effect of the cross-correlation parameter p^ on the exceedance density func- 
tion is illustrated in figure 7. It is this parameter which is neglected in the exceedance 
theories presented in references 2, 10, and 11. For the present choice of parameters the 
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Figure 6.- Exceedance density function showing effect of cross-correlation coefficient y. 


(For base case, °l, x - V = 0 1>x = U 2;X - 1; P - P' 0 = 0; Pi jX = M 2;X = 0; 
l 2c = >■) 


l lc 


exceedance density is a constant function throughout its domain 0 = 6 = 2ir. Reference to 
equations (71) to (78) reveals that this is indeed the case. As the cross-correlation 
parameter is increased (from p’ Q = 0) , the level of the density function is gradually 
reduced as illustrated in the figure. This same trend appears in the inset to figure 7, 
which shows the variation of the exceedances N + with p' Q . 

The case with p' Q = 0 has been discussed previously (fig. 4). It is of interest, 
however, to discuss the limiting case as p’ Q - 1.0. For p’ Q = 1.0, there is perfect cor- 
relation between x^ and Xg and between ^ and Xj. Therefore, with probability 1 



As a consequence of equations (80), the outward normal component V n of velocity 

Vn^XjXj+XgXg ( 81 ) 
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Figure 7.- Exceedance density function 
(For base case, cr^ - o 2x = 0 lx 

he ’ *•) 


showing effect of cross- correlation parameter 
‘Ax’ 1 ' P-r= 0; p 1iX = m 2jX = 0; x lc = 



vanishes with probability 1 and G + (0) = 0. This trend is indicated by the dashed line 
in the inset of figure 7. Therefore, the parameter p' Q , which is customarily nulled, may 
be Significant. 


CONCLUDING REMARKS 


A theory for combined random stresses has been developed for predicting the 
expected number of times per unit time that an n-dimensional random stress vector 
crosses a stress-interaction boundary (hyper surface). A joint statistical approach is 
employed and accounts for all cross-correlation effects between the stress components. 
An exact expression is obtained for the n-dimensional exceedance density function for an 
arbitrary stress boundary. 

The general n-dimensional theory is reduced for application to biaxial states of 
combined random stress. An elliptical stress-interaction boundary is assumed and the 
density function is expressed in a form convenient for calculating the exceedances by 
numerically evaluating a line integral. The characteristic behavior of the exceedance 
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density function for a circular stress boundary is discussed for some nominal parameter 
combinations. 

Of particular theoretical interest is the influence of a cross-correlation parameter 
p' Q which has usually t been neglected in previous exceedance theories. This parameter 
is representative of the measure of cross correlation between the stress and stress-rate 
components. For some combinations of parameters this cross- correlation effect was 
found to have a significant influence on the exceedances. It remains to be seen, however, 
whether the conditions under which this cross correlation becomes important are realized 
in practical combined, stress design applications. The improved exceedance theory pre- 
sented herein should provide a mathematical basis for further research in this area. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., December 28, 1973. 
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APPENDIX 


INTEGRALS OF QUADRATIC FORMS 


• For the convenience of the reader this appendix presents some useful properties 
relating to integration of quadratic forms. The notation follows very closely that of the 
main text and of reference 17. The same results can also be found elsewhere in the 
literature. 

It is desired to evaluate the integral 


p+OO p + OO „ + oo 

I = ... *(U) 

— oo ^-00 ^ _ 0 O 


du 


(Al) 


where the vector U has the components ^u^, Ug, 
$(U) = exp 


u n) 


and 


*\ 


1 T - 1 
U A C U 

2 




dU = dujdug. . .du n 


(A2) 


It will be assumed herein that C is a constant symmetric positive-definite matrix. 
There then exists a nonsingular n by n matrix Q such that 


C =Q T Q 


(A3) 



(A4) 


Now make a transformation of coordinates 


W - 



(A5) 


and denote the components of W as {wj, w 2 , . . ., w n }. By substituting equation (A4) 
into equation (Al), it follows that 


n+ QO n+°° 

p + °° 

2 2 

0 

M J ■ ' 

. \ exp 

Wi+w 2 + . . 

• + w n 

v — oo — OO 

j - oo 




J dwjdwg. 


,dw n 


(A6) 
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where J is the Jacobian determinant of the transformation (eq. (A5)). From the well- 
known result of reference 17, 


f” exp^r, 2 /^, , ^ 


(A7) 


and by using equation (A4), it follows that 


I = 



(A8) 


Partition U as 


U 


U = < 


U 





(b) 


J 


in which 


(A9) 


and 


r ~\ 
U 1 

u 2 






fu 


k+l) 


u k+2 


(AlOa) 


U 


(AlOb) 
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Consider the integral 

r>+°o /->+0O P+°° 

■M-U." -J.„ w)du (b) 


(All) 


where 



du k + l du k+2- • 


.diin 


The matrix C and the inverse 



are partitioned as follows: 


(A12) 


C = 


M 'T 

I 

t t :s 



(A13a) 


(A13b) 


The various submatrices appearing in equations (A13) have the following orders: 

M and Q - * k by k 

T and R k by (n - k) 

S and P"1 (n - k) by (n - k) 

The submatrices can be shown to satisfy the following relations (see, for example, 
ref. 17): 


pr t = -t t m _1 

(A14a) 

- 1 - 1 T 

Q = M + RPR 

(A 14b) 

P = S - t t m -1 t 

(A14c) 

M'MM ' 1 

(A14d) 
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With some manipulation, equation (All) can be written as 


1(b) = exp 


-Ka>M%) 

J 


p + oo p+OO p + 0O 

\ \ • • • \ exp 

‘■'-oo^'-oo -oo 


±lu (b) - V J p " 1 ( D (b) - v 


dU 


(b) 


(A15) 


in which 


V = -PR a U 


(a) 


(A16) 


Applying equation (A8) gives 


n-k 


I( b) = (2*) 


l£l 

M 


exp 


\ U fa) M_lu (a) 


(A17) 
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